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ABSTRACT 

To  maintain  the  continued  airworthiness  of  military  aircraft  it  is  essential  that  the 
fatigue  behaviour  of  components  subjected  to  complex  multi-axial  stress  conditions  be 
both  understood  and  predicted.  This  topic  is  extremely  complex.  Numerous  fatigue 
failure  criteria  ranging  from  the  purely  empirical  to  the  theoretical  have  been 
proposed.  These  criteria  rely  on  the  estimation  of  stress  and  strain  at  fatigue  critical 
locations.  This  interim  report  focused  on  possible  approaches  which  may  be  applicable 
to  both  low-  and  high-cycle  fatigue  regimes.  It  discusses  the  relative  advantages  of  the 
Neuber  and  the  Glinka  methods  for  calculating  localised  notch  strains  as  compared  to 
the  results  from  finite  element  analysis.  These  former  techniques  are  at  the  core  of 
several  sequence  accountable  crack  initiation  prediction  models,  some  of  which  are 
used  in  the  life  assessment  of  RAAF  aircraft.  Thus  the  accuracy  of  these  techniques 
directly  impact  the  estimated  fatigue  life  of  these  aircraft. 
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Stress  and  Strain  Estimation  at  Notches  in 
Aircraft  Structures 


Executive  Summary 

To  maintain  the  continued  airworthiness  of  military  aircraft  it  is  essential  that  the 
fatigue  behaviour  of  components  subjected  to  complex  multi-axial  stress  conditions  be 
both  understood  and  predicted.  Numerous  fatigue  failure  criteria  ranging  from  the 
purely  empirical  to  the  theoretical  have  been  proposed.  These  rely  on  the  accurate 
estimation  of  the  stress  and  strain  state  at  fatigue  critical  locations,  normally  associated 
with  stress  concentrations  in  the  form  of  notches.  This  interim  report  examines  the 
relative  advantages  of  the  Neuber  and  the  Glinka  methods  for  calculating  localised 
notch  strains.  These  former  techniques  are  at  the  core  of  several  sequence  accountable 
crack  initiation  prediction  models,  some  of  which  are  used  in  the  life  assessment  of 
RAAF  aircraft.  Thus  the  accuracy  of  these  techniques  directly  impact  the  estimated 
fatigue  life  of  these  aircraft.  The  accuracy  of  the  Neuber  and  Glinka  methods  was 
assessed  by  comparison  with  the  results  of  detailed  finite  element  analysis.  This  analysis 
confirmed  that  the  appropriateness  of  these  techniques  was  dependant  upon  the  stress 
state  (plane  stress  or  plane  strain)  at  the  notch  root. 
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Young's  modulus 
remote  strain 
normal  acceleration 

theoretical  elastic  stress  concentration  factor 
elastic-plastic  strain  concentration  factor 
elastic-plastic  stress  concentration  factor 
Ratio  of  minimum/  maximum  stress 
remote  stress 

Number  of  cycles  to  failure 
strain  energy  density 
strain 

elastic  strain 

strain  calculated  using  an  elastic-plastic  analysis  (total  strain) 
plastic  strain  contribution 

local  strain  components  in  plane  stress  conditions;  Cartesian  co-ordinates 

local  strain  components  in  plane  strain  conditions;  Cartesian  co-ordinates 

strain  at  the  yield  strength 

generalised  Poisson's  coefficient 

Poisson's  coefficient 

stress 

stress  calculated  using  an  elastic-plastic  analysis  (total  stress) 
local  stress  components  in  plane  stress  conditions;  Cartesian  co-ordinates 
local  stress  components  in  plane  strain  conditions;  Cartesian  co-ordinates 
yield  strength 
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1.  Introduction 


Aircraft  structural  components  are  generally  subjected  to  complex  loading  spectra.  The 
alternating  loads  tend  to  initiate  fatigue  cracks  at  notches  and  other  regions  of  high 
stress.  Historically  the  field  of  fatigue  has  been  classified  into  a  number  of  specific 
areas;  viz:  high-cycle  and  low-cycle  fatigue;  fatigue  of  notched  and  un-notched 
members;  the  initiation  and  propagation  of  cracks  and  fatigue  life  extension 
techniques.  Initially  the  focus  was  on  developing  simple  design  guidelines.  In  [1]  it  was 
remarked  that  in  the  case  of  high  cycle  fatigue,  the  fatigue  damage  was  associated  with 
localised  plasticity.  This  led  to  the  use  of  a  variety  of  approaches  based  on  the  concept 
of  an  equivalent  von  Mises  strain. 

In  this  interim  report  we  begin  by  briefly  reviewing  the  existing  strain  based  and 
energy  based  failure  criteria  theories.  Attention  is  focused  on  approaches  which  may 
be  applicable  to  both  low-  and  high-cycle  fatigue  regimes.  The  relative  advantages  of 
the  most  common  of  these  Neuber  [2]  and  the  Glinka  approaches  [3,  4]  for  calculating 
localised  notch  stresses  and  strains  are  discussed.  These  techniques  are  at  the  core  of 
several  sequence-accountable  crack  initiation  prediction  models,  some  of  which  are 
used  in  the  life  assessment  of  RAAF  aircraft. 

In  a  recent  review  [5]  of  elastic-plastic  notch  root  stress-strain  estimation  methods,  the 
following  general  conclusions  as  to  the  accuracy  of  Neuber  and  Glinka  estimations 
were  drawn: 

•  The  relative  error  of  the  Neuber's  rule  is  in  general  larger  than  that  of  Glinka. 
However,  Glinka  normally  over-estimates  where  Neuber  under-estimates  the  notch 
root  strain.  Neuber  will  therefore  be  generally  conservative  while  Glinka  is  non¬ 
conservative. 

•  The  closer  the  loading  situation  is  to  the  elastic  case,  the  smaller  is  the  relative  error 
between  the  two  methods. 

•  Neuber's  method  generally  gives  the  best  prediction  under  predominantly  plane 
stress  state,  whilst  Glinka's  method  is  best  for  nearly  plane  strain  conditions. 

These  generalizations  are  investigated  in  this  report  by  comparing  the  estimations 
based  on  the  Neuber  and  Glinka  methods  with  die  results  of  a  detailed  finite  element 
analysis.  For  this  purpose  a  rectangular  plate  containing  a  circular  hole  was 
investigated.  Both  plane  stress  and  plane  strain  conditions  were  investigated  for  two 
levels  of  monotonic  loading. 

A  number  of  recommendations  are  then  made  to  address  some  of  the  shortcomings 
inherent  in  the  current  crack  initiation  lifing  methodology. 
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2.  Common  Fatigue  Failure  Criteria 


2.1  HIGH  CYCLE  FATIGUE:  Stress  Based  Criteria 

Stress  based  methods  are  usually  associated  with  high  cycle  fatigue.  Whilst  there  are  a 
large  variety  of  failure  criteria  used  many  of  these  criteria  can  be  written  in  the  form; 
viz: 


f(I,,I2)  =  q(Nf)  (1) 

where  Ii  =  oi  +  a2+  03  and  I2  =  ai  a2+  c2  03+  01  03 ,  are  the  first  and  second  invariants  of 
the  stress  tensor.  Here  ai ,  o2  and  03  are  the  principal  stresses,  01  >  a2  >  03  and  Nf  is 
number  of  cycles  to  failure.  One  commonly  used  form  for  q  is 

q(Nf)  =  K  N“  +  C  (2) 

where  K  is  a  constant,  a  <  0  and  C  >  0.  If  Nf  00  then  eqn.  (1)  becomes 

f(Ii,I  i)  =  C  (3) 

which  resembles  the  yield  criteria  for  isotropic  metals.  One  of  the  most  familiar,  and 
widely  used,  yield  criteria  is: 


1/2  Sy  Sy  =C0  (4) 

where  Sy  is  the  deviatoric  stress  tensor,  Co  =  3oy  and  ay  is  the  uniaxial  tensile  yield 
stress.  Many  authors  have  adopted  this  analogy  and  some  of  the  commonly  used  failure 
(yield)  criteria  are  given  in  [6,  7,  8,  9].  Suitable  adjustment  of  these  criteria  can  reduce 
the  majority  to  either  the  Tresca  or  the  von  Mises'  yield  condition.  The  fact  that  stress 
based  criteria  are  more  applicable  to  high  cycle  fatigue  is  discussed  in  [10]. 


2.2  LOW  CYCLE  FATIGUE:  Strain-based  criteria 

In  contrast  strain  based  methods  are  usually  associated  with  low  cycle  fatigue.  In  this 
post  initial  yield  regime  the  Manson-Coffin  relationship  [1 1, 12]: 


As/2  =  Ase  /2  +  Asp  /2  =  af  /E  (Nf)b  +  sf'  (Nf) c  (5) 

is  widely  used.  Here  A  refers  to  its  ’range1;  and  cxf  /E  and  s  f '  are  strain  amplitudes 
corresponding  to  the  elastic  and  the  plastic  intercept  for  one  cycle  respectively;  E  is 
Young's  modulus,  and  b  and  c  are  experimental  constants. 
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To  extend  this  law  to  the  multi-axial  stress  states  requires  the  definition  of  an 
equivalent  strain.  This  extension  is  often  written  in  the  form 

Aseq/  2  =  Of  /  E  (Nf  )b  +  s  f ' ‘  (Nf) c  (6) 

A  number  of  expressions  for  the  equivalent  strains  have  been  used;  viz:  the  octahedral 
shear  strain  [10,  13],  the  maximum  shear  strain  [10,  14],  the  maximum  normal  strain 
[10, 15],  the  von  Mises'  equivalent  strain  [16],  modified  von  Mises'  [14, 17, 18]  and  the 
maximum  total  strain  [19]. 

One  expression  commonly  used  as  an  equivalent  strain  measure  is: 

6eq  =  a((ei-S2)2+(E3-£2)2+  (Sl-Es)2)1^  (7) 

When  a=  1/3  we  find  that  e  eq  =y0ct,  the  octahedral  shearing  strain.  When  a=  ^2/ 2(1+  o) 
then  Eeq  coincides  with  the  expression  for  the  von  Mises  equivalent  strain. 

3.  Notch  Root  Stresses  and  Strains 

The  strain-life  fatigue  method  requires  that  the  notch  root  stresses  and  strains  be 
known.  These  quantities  can  be  determined  in  several  ways,  viz:  via  direct  strain 
(mirco)  measurements,  using  finite  element  analysis  or  by  using  approximate  methods 
that  relate  local  stresses  and  strains  to  their  remote  values. 

The  theoretical  stress  concentration  factor,  Kt/  is  often  used  to  relate  tire  nominal  (or  far 
field)  stresses  S,  or  strains,  e,  to  the  local  values,  a  and  s.  Upon  yielding,  Hooke's  law 
cannot  be  used  to  relate  the  local  stress,  a,  to  the  local  strain,  e,  and  the  local  values  are 
no  longer  related  to  the  nominal  values  by  Kt  Instead,  the  local  stresses  and  strains  are 
related  to  the  remote  values  by  their  respective  stress  and  strain  concentration  factors 
Ka  and  Ke ,  viz: 


a  =  Ko  S  and  s  =  KE  e  (8) 

In  an  attempt  to  compute  the  local  (notch)  stress  and  strains  the  strain  life  approach 
frequently  makes  use  of  Neuber's  rule  [2].  In  this  approach  the  theoretical  stress 
concentration  is  taken  as  the  geometric  mean  of  the  stress  and  strain  concentration 
factors.  It  is  widely  assumed  that  this  relationship  holds  true  for  most  notch  geometries 
(see  Potter  in  [20]).  Problems  associated  with  this  and  other  related  approximations  are: 

1.  The  notch  root  stress-strain  response  may  not  always  be  in  phase  with  the  global 
load  [21]; 

2.  The  law  does  not  account  for  time  dependent  processes  such  as  creep  and  stress 
relaxation;  and 
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3.  The  law  does  not  account  for  cyclic  stress  relaxation. 

A  number  of  variants  of  Neuber's  rule  have  been  used  to  relate  the  remote  stresses  and 
strains  to  local  values.  The  basic  Neuber  approach  and  a  commonly  used  variant  is 
outlined  below: 

The  initial  Neuber  hypothesis  [2]  can  be  expressed  in  the  form: 

K,2  =  Ka  Ke  (9) 


or  alternatively 


S  e  Kt  2  =  a  e  (10) 

When  the  remote  stresses  and  strains  are  above  yield  Seeger  and  Heuler  [22] 
proposed  an  extension  of  Neuber's  rule,  viz: 

SeKp2  =  g  e  (11) 

where 

Kp  =  Stress  S  at  the  onset  of  general  yielding  /  Stress  S  at  which  the  notch  first  yields 
and  S_=  S  Kt  /  Kp  (12) 

Here  Sand_e  must  lie  on  the  cyclic  stress-strain  curve. 

The  formulation  of  equation  (10)  is  widely  accepted  [20]  and  is  used  in  many 
current  sequence  accountable  crack  initiation  prediction  models,  see  [23,  24,  25, 
26,  27].  Currently  the  F/  A-18  structure  is  monitored  using  a  derivative  of  one  of 
these  programs. 

In  the  authors  opinion  these  formulations  have  a  number  of  shortcomings;  viz: 

i)  For  cyclic  loading  it  is  common  to  modify  Neuber's  hypothesis  by  replacing  Kt  by  a 
fatigue  notch  reduction  factor  K/  see  Topper  et  al.  [28]. 

The  work  of  Glinka  [3,  4]  has  shown  that  Neuber’s  hypothesis  can  incorrectly  estimate 
the  inelastic  strain.  The  need  to  replace  Kt  by  Kf  etc,  can  be  overcome  by  using  the 
Glinka  hypothesis  to  calculate  As.  According  to  this  hypothesis  the  strain  energy  density 
at  the  notch  in  a  fully  plastic  analysis  (small  scale  plasticity)  equates  to  the  strain  energy 
obtained  via  a  purely  elastic  analysis;  viz: 


1/2  Kt  Se  =  dsy  =  la  ds  (13) 
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Once  a  valid  stress  strain  relationship  is  known  the  notch  stresses  and  strains  can  be 
readily  evaluated. 

ii)  For  plane  stress  or  plane  strain  problems,  one  approach  is  to  use  simple  power  laws 
to  relate  a  to  s,  or  for  saturated  fatigue  loops  Act  to  As.  These  are  normally  formulated 
assuming  the  material  exhibits  Masing’s  hypothesis  [31],  which  states  that  following 
initial  yielding,  the  shape  of  the  hysteresis  curve  can  be  approximated  by  assuming  it  is 
similar  to  the  initial  loading  curve  magnified  by  a  factor  of  two.  This  assumption  is  not 
appropriate  for  all  materials. 

Alternatively,  the  stress-strain  relationships  can  be  “digitised”  as  an  input  into  the  lifing 
models  [see  23, 24, 25, 26, 27]. 

iii)  The  relations  presented  above,  and  as  used  to  life  some  military  aircraft  (eg  the  F/A- 
18),  do  not  allow  for  the  local  level  of  constraint/triaxiality.  One  method  for  including 
constraint  effects,  which  modifies  the  expression  for  the  constant  K  in  the  fatigue  law,  is 
outlined  in  [29]. 

For  complex  3D  structural  problems  the  use  of  a  validated  constitutive  law  may  be 
necessary.  In  this  case  when  calculating  the  Ag  and  As  ranges  the  use  of  traditional 
yield  surface  plasticity  should  be  avoided  as  this  tends  to  produce  “boxy”  stress  strain 
loops  and  gives  poor  estimates  of  the  Ag  and  As  ranges.  Stiffener  runout  Number  2  in 
the  FI  1 1  wing  pivot  fitting  is  a  good  example  of  this  problem.  Here  when  subjected  to  g 
loads  going  from  0 g  to  7.3g  to  Og  simple  classical  incremental  plasticity  gave  a  tensile 
stress  range  of  approximately  2,000  MPa  whilst  more  exact  analysis,  using  a 
constitutive  law,  gave  a  stress  range  of  only  1,200  MPa  [30], 

vi)  Care  should  be  taken  to  ensure  these  relations  adequately  allow  for  the  load  history 
effects,  i.e.  load  interaction,  stress  relaxation,  creep,  and  overload-underloads. 

Here  the  stiffener  runout  Number  2  in  the  Fill  wing  pivot  fitting  is  (again)  a  good 
example  of  the  need  to  correctly  follow  load  history.  When  subjected  to  g  loads  going 
from  0 g  to  7.3 g  to  0 g  simple  classical  approaches  gave  inspection  interval  of  less  than 
500  hours.  However,  allowing  for  the  load  history  enabled  the  inspection  interval  to  be 
increased  to  almost  1,500  hours  [30]. 

4.  Comparison  of  Neuber's,  Glinka's  and  Finite 

Element  Results 


As  we  have  previously  seen  the  accurate  prediction  of  the  stresses  and  the  strains  in  a 
structure  are  important  in  fracture  mechanics  and  fatigue  life  prediction.  In  this  section 
of  the  report  we  will  evaluate  the  relative  ability  of  Neuber's  rule  and  Glinka's  method 
to  calculate  the  local  stresses  and  strains  in  a  plate  with  a  circular  hole.  To  this  end  the 
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values  predicted  via  Neuber's  rule  and  Glinka's  method  were  compared  against  finite 
element  method  (FEM)  results  obtained  using  classical  incremental  plasticity. 

As  an  example  let  us  consider  a  22mm  wide  x  64mm  long  rectangular  plate,  1mm  thick 
with  a  8mm  diameter  circular  hole  located  at  the  centre  of  the  plate,  see  Figure  1.  The 
plate  is  subjected  to  a  remote  uniform  stress.  The  load  spectra  applied  to  the  plate 
involved  monotonic  loading  from  zero  to  a  peak  value  of  a,  the  load  was  then  fully 
reversed  until  a  remote  stress  of  -cr  was  reached,  at  which  stage  the  load  was  again 
reversed  until  a  remote  stress  of  crwas  reached,  ie.  the  cycle  was  (0,  cr,  -cr,  cr),  R  =  -1 . 


o 


Figure  1:  Schematic  diagram  of  a  rectangular  plate  with  a  circular  hole 

This  analysis  was  repeated  with  two  different  levels  of  cr ;  viz:  cr  =  150  MPa  and  225 
MPa.  These  two  values  of  the  remote  stress  were  selected  so  that  the  nominal  section 
would  be  at  different  levels  of  stress  which  both  approach  the  yield  strength  of  the 
material.  The  load  was  applied  in  increments  so  the  cyclic  behaviour  of  the  plate  could 
be  observed  and  the  plate  was  considered  to  be  in  a  state  of  either  plane  stress  or  plane 
strain.  In  this  analysis  the  yield  strength,  crys,  of  the  idealised  material  was  taken  as  400 
MPa,  with  a  Young's  Modulus,  E,  of  70,000  MPa  and  a  post  yield  gradient  of  2,003 
MPa  (elastic  -  linear  strain  hardening).  The  stresses  and  plastic  strains  at  the  notch  were 
the  calculated  using  the  FEM,  Neuber's  rule  and  Glinka's  method  at  node  2  as  depicted 
in  Figure  1. 

The  plate  was  modelled  using  8-noded  plate  elements.  A  denser  mesh  was  used 
around  the  hole  to  model  the  stresses  and  strains  to  a  greater  accuracy  since  steeper 
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stress  gradients  are  present  and  the  material  yields  at  and  in  the  vicinity  of  the  notch 
tip.  Because  of  symmetry  only  a  quarter  of  the  plate  was  modelled,  a  typical  mesh  is 
shown  in  Figure  2. 


Figure  2:  A  typical  mesh  of  the  modelled  plate 


4.1  Neuber's  Rule:  -  Plane  Stress 

For  plane  stress  there  is  only  a  uniaxial  stress,  cry  at  node  2  since  at  this  point  crx-  cz-  0. 
The  monotonic  stress-strain  curve  used  for  the  material  is  given  by  equations  (14)  and 

(15). 


=  Ee  if  cr  <  400 MPa 

(14) 

=  2003^  +  388.55  if  cr>400 

(15) 

The  cyclic  hysteresis  stress-strain  curve  used  is  given  by  equations  (16)  and  (17).  This 
example  assumed  that  the  stress  -  strain  hysteresis  curve  is  approximately  twice  the 
value  of  the  cyclic  stress  -  strain  curve  (Masing's  Hypothesis  [31]).  This  assumption  has 
a  sound  basis  since  the  cyclic  stress  -  strain  curve  is  the  locus  of  the  tips  of  the 
hysteresis  loops.  When  the  material  has  been  taken  past  the  yield  point,  the  new  origin, 
in  stress  strain  space,  is  located  at  the  last  value  of  the  stress  and  strain  at  which  a  load 
reversal  occurred. 
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-  Es  if  a  <  800 MPa 


(16) 


=  2003f  +  777.1 1  if  cr  >800 MPa  (17) 

Using  Neuber's  approach  and  the  stress-strain  curve  for  monotonic  loading  allows  the 
values  of  stress,  crf  and  the  plastic  strain,  et  to  be  evaluated.  Similarly  using  Neuber's 
approach  and  the  cyclic  hysteresis  stress-strain  curve  for  cyclic  loading  allows  the 
values  of  stress,  crt  and  the  strain,  &  to  be  evaluated  at  any  point  in  the  load  cycle.  Here 
the  subscript  "t"  indicates  that  the  values  were  calculated  using  an  elastic-plastic  stress 
analysis.  For  monotonic  loading  Neuber's  approach  gives: 


2003£,2  +388.55^  - 


=  0 


if 


a  >  400 MPa 


(18) 


Solving  this  quadratic  equation  allows  st  to  be  evaluated. 

A  similar  substitution,  using  the  cyclic  hysteresis  stress  strain  curve,  can  be  made  for 
cyclic  loading  to  solve  for  et.  For  the  cyclic  hysteresis  stress  strain  curve  it  was 
assumed,  for  calculation  purposes,  that  the  origin  was  placed  at  the  last  stress  and 
strain  value  before  the  load  was  reversed  and  the  positive  axes  of  stress  and  strain  are 
in  the  direction  of  the  reversed  load. 

4.2  Glinka's  Method  -  Plane  Stress 

When  using  Glinka's  method  together  with  the  stress  strain  curve  for  monotonic 
loading  we  obtain  a  similar  quadratic  equation  allowing  the  stress,  crf  and  the  strain,  et 
to  be  calculated;  viz: 


0  =  2003£’,2  +  £,(388.55  +  cr  -  2003 e  )  + 


if  cr  >400 MPa  (=  crys) 


(19) 


For  cyclic  loading  use  of  the  cyclic  hysteresis  stress  strain  curve,  equation  (17)  again 
gives  a  quadratic  equation  allowing  the  stress,  ot  and  the  strain,  £>to  be  evaluated. 
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4.3  Neuber's  Rule  and  Glinka's  Method  in  Plane  Strain 

For  plane  strain,  an  analogous  analysis  can  be  carried  out  as  per  that  outlined  above  for 
the  plane  stress.  However,  in  this  case  a  biaxial  stress  state  exists  at  the  notch  tip 
(because  oy'  &  0,  oz'  #  0).  To  allow  for  this  Dowling  et  al.  [32]  derived  the  following 
relationships  which  allow  for  the  translation  of  the  uniaxial  stress  strain  curve  into  a 
biaxial  "plane  strain"  stress  stain  relation;  viz: 


(  \ 


(  ,  ^ 

I  1  I 


(v  +  E£p  llo  )  o 

where  n  =  - - — r- ,  £  =  —  +  =  st  +  sv 

(l  +  Efp/rr)  E 


(20) 


(21) 

(22) 


Thus,  under  plane  strain  conditions,  the  same  set  of  equations  using  Neuber's  rule  and 
Glinka's  method  can  be  used  except  that  all  the  calculations  are  based  on  the  "plane 
strain"  stress-strain  curve  cry 


4.4  Discussion 


The  accuracy  of  the  FEM  model  was  first  checked  against  the  elastic  theoretical 
solution  at  the  notch  tip  when  the  notch  tip  was  still  within  the  linear  elastic  range 
during  monotonic  loading  of  the  plate.  The  theoretical  solution  at  the  notch  tip  when 
the  material  is  linear  elastic  is  given  by. 


and 


|  y  ‘  for  plane  stress 

Gy  —  Kt(T„om 

)  <7^=0  for  plane  strain 

o'-  vK.o„nm 

k  z  r  nom 


(23) 


(24) 
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In  this  case  Peterson  [33]  gave  a  Kt  of  3.658.  The  computed  values  for  all  the  stress 
components  are  shown  in  Table  3.  In  each  case  the  numerical  results  were  within  2%  of 
the  theoretical  value. 

Table  3:  Percentage  error  in  the  FEM  model 


(For  node  2) 


FEM 

THEORETICAL 

FEM 

THEORETICAL 

PLANE 
STRAIN 
(%  error) 

PLANE 

STRAIN 

PLANE 
STRESS 
(%  error) 

PLANE  STRESS 

CTx(MPa) 

0.01 

0 

0.02 

0 

ay(MPa) 

1 

400.55 

WKmlM 

406.04 

oz(MPa) 

HIRBhH 

128.17 

0 

0 

4.4.1  Comparison  of  Neuber's  Rule  and  Glinka's  method  to  the  FEM  model 

Having  confirmed  the  accuracy  of  the  model  the  computed  hysteresis  loops,  for  both 
plane  stress  and  plane  strain,  and  for  both  load  cases  are  shown  in  Figures  3  to  6.  From 
these  figures  it  can  be  observed  that  all  three  methods  predicted  a  similar  cyclic 
behaviour  at  the  notch  tip  but  with  different  magnitudes  for  the  strains.  All  models 
displayed  elastic  unloading  and  symmetry  in  the  cyclic  stress-strain  curve.  Each 
method  also  predicted  "boxy"  stress-strain  curves,  as  would  be  expected  given  that  the 
theory  for  each  method  was  based  on  incremental  plasticity  and  an  idealised  material. 
In  reality,  "curved"  stress-strain  curves  should  be  expected. 

The  percentage  error  of  Neuber's  rule  and  Glinka's  method,  as  compared  with  the  FEM 
results,  in  predicting  the  maximum  stresses  and  strains  in  monotonic  loading  is  shown 
in  Table  4. 


Table  4:  The  percentage  error  in  the  maximum  stresses  and  strains  under  monotonic 

loading  using  Glinka's  method  and  Neuber's  rule 


Load 

Type  of 
Loading 

FEM  model 

Glinka's  Method 
(%  error) 

Neuber's  rule 
(%  error) 

MPa 

e 

a  (MPa) 

s(%) 

G(%) 

6<%) 

a  (%) 

150 

Plane 

Stress 

0.00867 

406.57 

0.00821 

(-5.22%) 

405.01 

(-0.38%) 

0.01050 

(21.1%) 

409.59 

(0.74%) 

225 

Plane 

Stress 

0.01519 

424.64 

0.01474 

(-2.94%) 

418.09 

(-1.54%) 

0.02233 

(46.9%) 

433.29 

(2.04%) 

150 

Plane 

Strain 

0.00721 

456.13 

0.00761 

(5.54%) 

460.91 

(1.05%) 

0.00918 

(27.3%) 

466.71 

(2.32%) 

225 

Plane 

Strain 

0.01247 

476.00 

0.01330 

(6.67%) 

479.43 

(0.72%) 

0.01945 

(55.9%) 

496.55 

(4.31%) 
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4.4.2  Remarks 

For  the  particular  problem  studied  this  analysis  has  shown  that,  for  the  material  stress- 

strain  law,  geometry  and  loading  considered,  Glinka's  method  estimates  the  stresses 

and  strains  to  a  greater  accuracy  than  Neuber's  rule.  The  following  general 

observations  can  be  made. 

1.  Both  Neuber's  rule  and  Glinka's  method  calculated  the  stresses,  oi,  to  within  5%  of 
the  FEM  results. 

2.  Neuber's  rule  significantly  overestimated  the  strain,  St,  for  both  load  cases  and  both 
load  types. 

3.  As  the  nominal  stress  was  increased,  the  errors  in  the  predictions  made  by 
Neuber's  rule  increased  significantly. 

4.  For  plane  stress  Glinka's  method  underestimates  the  value  of  the  strain.  However, 
it  is  within  approximately  5%  of  the  calculated  values.  For  plane  strain  Glinka's 
method  overestimates  the  value  of  strain  and  stress  but  is  within  approximately  5% 
of  the  calculated  value. 

5.  As  the  nominal  stress  in  the  plate  was  increased  (but  was  still  below  the  yield 
strength  of  the  material)  Glinka's  method  predicted  the  stresses  and  strain  to 
approximately  the  same  degree  of  accuracy  as  the  FEM. 

These  findings  are  generally  in  agreement  with  those  of  Shin  [5]. 
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Stress  -  Strain  Curve  for  150  MPa  loading  in  Plane  Stress 

1  1 

1  1 

1  1 

1  ! 

1  1 

1  1 

/ 

l  SJ.OOOOj 

Figure  3:  Stress  Strain  Curves  for  155  MPa  (Plane  Stress) 

Stress  -  Strain  Curve  for  225  MPa  loading  in  Plane  Stress 
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Figure  4:  Stress  Strain  Curves  for  22 5  MPa  (Plane  Stress) 
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ure  5:  Stress  Strain  Curves  for  155  MPa  (Plane  Strain) 


Stress  -  Strain  Curve  for  225MPa  loading  in  Plane  Strain 
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Figure  6:  Stress  Strain  Curve  for  250  MPa  (Plane  Strain) 


DSTO-GD-0196 


COE-SM-98-Ol 


5.  Conclusion 

This  report  has  presented  a  brief  summary  of  some  common  relationships  used  to 
estimate  notch  stress  and  strain  for  both  low-  and  high-cycle  fatigue.  Particular 
attention  has  been  paid  to  the  relative  advantages  of  the  Neuber  and  the  Glinka 
approaches. 

Current  lifing  methodology  is  frequently  based  on  the  an  early  version  of  the  strain  life 
approach  and  has  a  number  of  shortcomings  which  should  be  addressed;  viz: 

1.  Neuber' s  hypothesis  is  frequently  used  to  compute  the  strain  range.  When  using 
this  approach  Kt  should  be  replaced  by  the  fatigue  notch  reduction  factor  K/ ,  see 
Topper  et  al  [28],  However,  it  is  desirable  that  both  Neuber's  hypothesis  and 
Glinka's  hypothesis  should  be  available  for  use  in  any  lifing  analysis  package.  This 
would  require  the  user  to  have  a  prior  knowledge  of  the  relative  accuracy  of  each 
method  for  the  particular  material,  geometry  and  loading  spectrum  under 
consideration. 

2.  Some  analyses  use  power  laws  to  compute  Aa  and  Ac.  Some  relationships  used  are 
only  accurate  for  Masing  type  materials.  The  use  of  generalisations  of  this 
approximation  for  non  Masing  materials  should  be  investigated. 

3.  The  strain  life  notch  relations,  as  currently  used,  do  not  allow  for  the  local  level  of 
constraint/  triaxiality.  This  can  have  a  very  significant  effect.  Methods  for  including 
constraint  effects  should  be  adopted. 

4.  Care  should  be  taken  to  ensure  these  relations  adequately  allow  for  the  load  history 
effects,  i.e.  load  interaction,  stress  relaxation,  creep,  and  overloa d-underloads . 
General  ways  of  accounting  for  these  effects  should  be  included  in  the  lifing 
formulae. 

In  comparison  with  the  results  of  a  finite  element  analysis  of  a  plate  with  a  central  hole: 

5.  Neuber's  rule  significantly  overestimated  the  notch  strain  for  both  plane  stress  and 
plane  strain  for  the  two  load  cases  considered. 

6.  As  the  nominal  stress  was  increased,  the  errors  in  the  predictions  made  by 
Neuber's  rule  increase  significantly. 

7.  For  plane  stress  Glinka's  method  underestimated  the  value  of  the  strain.  However, 
it  is  within  approximately  5%  of  the  calculated  values.  For  plane  strain  Glinka's 
method  overestimated  the  value  of  strain  and  stress  but  is  within  approximately  5% 
of  the  calculated  value. 
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8.  As  the  nominal  stress  in  the  plate  was  increased  (but  was  still  below  the  yield 
strength  of  the  material)  Glinka's  method  predicted  the  stresses  and  strain  to 
approximately  the  same  degree  of  accuracy  as  the  finite  element  method. 


5.1  Summary 

The  accuracy  of  the  current  analysis  methodology  generally  accepted  to  compute  the 
fatigue  (initiation)  life  may  have  the  potential  to  be  dramatically  improved.  Before 
conducting  a  lifing  analysis,  the  relative  accuracy  of  the  Neuber  or  Glinka  hypotheses 
for  the  particular  geometry  under  consideration  should  be  assessed,  and  the 
appropriate  hypothesis  chosen.  Further  attention  should  be  given  to  accounting  for  the 
level  of  local  constraint,  and  on  the  use  of  methodologies  capable  of  accounting  for  the 
effect  of  overloads  and  stress  relaxation. 

6.  Acknowledgments 

This  work  was  conducted  for  DSTO's  Airframes  and  Engines  Division  under  the 
auspices  of  the  DSTO  Centre  of  Expertise  -  Structural  Mechanics.  The  authors  wish  to 
acknowledge  the  support  and  useful  discussions  held  with  the  COE  liaison  officers  Drs 
A.  A.  Baker  and  L.F.R.  Rose,  and  C.  H.  Wang  of  AED. 

7.  References 


1.  Guest,  J.  J.,  "The  problem  of  combined  stress".  Engineering,  London,  155  (Jan.-Apr. 
1943)  21-3, 101-2, 281-2, 303-4. 

2.  Neuber  H.,  "Theory  of  stress  concentration  for  shear  strained  prismatic  bodies  with 
arbitrary  nonlinear  stress  strain  law".  Journal  of  Applied  Mechanics  28,  544-551, 
1969. 

3.  Glinka  G.,  "Energy  density  approach  to  calculation  of  inelastic  stress-strain  near 
notches  and  cracks".  Engineering  Fracture  Mechanics,  22, 3, 485-508, 1985. 

4.  Glinka  G.,  "Calculation  of  inelastic  notch-tip  strain-stress  histories  under  cyclic 
loading".  Engineering  Fracture  Mechanics,  22, 5, 839-854, 1985. 

5.  Shin,  C.S.,  "Fatigue  crack  growth  from  stress  concentrations  and  fatigue  life 
prediction  in  notched  components",  A.  Capinteri  (Editor),  Handbook  of  Fatigue 
Crack  Propagation  in  Metallic  Structures,  Elsevier  Science  B.V.,  1994. 

6.  Gough,  H.  J.  &  Pollard,  H.  V.,  "Strength  of  metals  under  combined  alternating 
stress",  Proc.  Inst.  Mech.  Engrs.,  131,  pp  3-54, 1935. 


15 


DSTO-GD-0196 


COE-SM-98-Ol 


7.  Stanfield,  G.,  Proc.  Inst.  Mech.  Engrs.,  131,  pp  59, 1935. 

8.  Findley,  W.  N.,  "Fatigue  of  metals  under  combinations  of  stress",  Trans.  ASME,  79, 
pp  1337-48, 1957. 

9.  McDiarmid,  D.  L.,  "A  new  analysis  of  fatigue  under  combined  bending  and 
twisting".  Aero.  /.,  Roy  Aero.  Soc.  London,  78  (1974)  325-9. 


10.  Sines,  G.  &  Ohgi,  G.,  "Fatigue  criteria  under  combined  stresses  or  strains",  /.  Engng. 
Mater.  &  Tech.  Trans.  ASME,  103,  pp  82-90, 1981 

11.  Tavemelli,  J.  F.  &  Coffin,  L.  F.  Jr,  "Experimental  support  for  generalized  equation 
predicting  low  cycle  fatigue",  ].  Basic  Engng.,  Trans.  ASME,  84,  pp  533-537, 1962. 

12.  Manson,  S.  S.,  J.  Basic  Engng.,  Trans.  ASME,  84,  pp  537-41, 1962. 

13.  Zamrik,  S.  Y.  &  Goto,  T.,  "The  use  of  octahedral  shear  strain  in  biaxial  low  cycle 
fatigue".  In  Materials  Technology  on  Inter-American  Approach.  ASME,  New  York, 
1968,  pp.  551-62. 

14.  Andrews,  J.  M.  H.  &  Ellison,  E.  G.,  "A  testing  rig  for  cycling  at  high  biaxial  strains", 
J.  Strain  Anal.,  8 ,  pp  168-75, 1973. 

15.  Libertiny,  G.  Z.,  "Short-life  fatigue  under  combined  stress",  J.  Strain  Anal,  2,  91-5, 
1967. 

16.  Pascoe,  K.  J.  &  DeVilliers,  J.  W.  R.,  "Low  cycle  fatigue  of  steels  under  biaxial 
straining",  J.  Strain  Anal.,  2,  117-26, 1967. 

17.  Havard,  D.  G.,  Williams,  D.  P.,  &  Topper,  T.  H„  "Biaxial  fatigue  of  mild  steel:  data 
synthesis  and  interpretation",  Ontaria  Hydro.  Res.  Quart.,  27,  pp  11-18, 1975. 

18.  Lefebvre,  D.  F.,  "Hydrostatic  effect  on  the  life  prediction  in  biaxial  low-cycle 
fatigue",  Proc.  2nd  Int.  Conf.  on  Multi-axial  Fatigue,  Sheffield  UK,  1985. 

19.  Zamrik,  S.  Y.  &  Frismuth,  R.  E.,  "The  effect  of  out  of  phase  biaxial  strain  cycling  on 
low-cycle  fatigue",  Exp.  Mech.  SESA,  13,  pp  204-8, 1973. 

20.  "Cyclic  stress-strain  behaviour:  Analysis,  experimentation,  and  failure  prediction", 
ASTM,  Special  Technical  Publication  519,  Philadelphia,  Pa.,  1972. 

21.  Paul,  J,  and  Molent,  L.,  "Applications  of  energy  density  theory  in  cyclic  plasticity",  J. 
Theoretical  and  Applied  Frac.  Mech.  10  (1988)  pp.43-48. 


16 


COE-SM-98-Ol 


DSTO-GD-0196 


22.  T.  Seeger  and  Huler  P.,  "Generalised  application  of  Neuber's  rule",  Journal  of  Test 
Eval.,  8, 199-204, 1980. 

23.  "Computer  program  documentation  for  C_CI89.",  McDonnell  Douglas  Aircraft,  St 
Louis,  USA,  1991. 

24.  "Introduction  to  fatigue  life  program  'CRACKSTART'.",  IFOSTP  Memorandum, 
Canadair,  Canada,  1990. 

25.  Porter,  P.G.  and  Liu,  A.F.,  "A  rapid  method  to  predict  fatigue  crack  initiation,  Vol  1 
&  2  -  Technical  Summary",  N ADC-8101 0-60,  Pa.  USA,  Feb  1983. 

26.  Ghidella,  J.R.  and  Madley,  W.B.,  "A  comprehensive  review  of  a  fatigue  analysis 
program  -  CI_89",  DSTO-TR  (in  press),  Melbourne,  Aust. 

27.  "CIJLIFE  user's  manual  for  the  prediction  of  crack  initiation  life",  MDC  A8873, 
McDonnell  Douglas  Aircraft,  St  Louis  USA,  Aug  1984. 

28.  Topper  T.  FI.,  Wetzel  and  Morrow  J.  D.,  "Neuber's  rule  applied  to  fatigue  of 
notched  specimens",  /.  Material  Science,  pp  200-209, 1969. 

29.  Golos,  K.  &  Ellyin,  F.,  "Cumulative  damage  and  fatigue  of  A-516  Gr.  70  pressure 
vessel  steel",  J.  Press.  Vessel  Tech.  Trans.  ASME,  vol  110,  pp  36-41,  1988. 

30.  R.  Jones,  L.  Molent,  J.  Paul,  T.  Saunders,  W.K.  Chiu,  "Development  of  a  composite 
repair  and  the  associated  inspection  intervals  for  the  F111C  stiffener  runout  region", 
FAA/NASA  International  Symposium  on  Advanced  Structural  Integrity  Methods  for 
Airframes  Durability  and  Damage  Tolerance,  NASA  Conference  Publication  3274,  Part 
1,  pp  339-351, 1994. 

31.  Masing,  G.,  "Internal  Stresses  and  Hardening  of  Brass",  Proc.  2nd  Int.  Congress  for 
Appl.  Mech.,  Zurich,  Sept  1926. 

32.  N.E.  Dowling,  W.R.  Brose  and  W.K.  Wilson,  "Notched  member  fatigue  life 
predictions  by  the  local  strain  approach".  In  Advances  in  Engineering,  vol.  6  (Edited 
by  R.M.  Wetzel),  p55-84.  Society  of  Automotive  Engineers,  SAE  Warrendale,  Pa. 
(1979). 

33.  R.E.  Peterson,  "Stress  Concentration  Design  Factors",  John  Wiley  &  Sons,  Inc,  New 
York,  p86  (1953). 


17 


DSTO-GD-0196 


COE-SM-98-Ol 


18 


DISTRIBUTION  LIST 


Stress  and  Strain  Estimation  at  Notches  in  Aircraft  Structures 

R.  Jones,  M.  Knopp,  J.  Price  and  L.  Molent 

AUSTRALIA 


DEFENCE  ORGANISATION 

Task  Sponsor 
CLSA 

S&T  Program 

Chief  Defence  Scientist 
FAS  Science  Policy 
AS  Science  Corporate  Management 
Director  General  Science  Policy  Development 
Counsellor  Defence  Science,  London  (Doc  Data  Sheet ) 

Counsellor  Defence  Science,  Washington  (Doc  Data  Sheet ) 

Scientific  Adviser  to  MRDC  Thailand  (Doc  Data  Sheet ) 

Director  General  Scientific  Advisers  and  Trials/Scientific  Adviser  Policy  and 
Command  (shared  copy) 

Navy  Scientific  Adviser  (Doc  Data  Sheet  and  distribution  list  only) 

Scientific  Adviser  -  Army  (Doc  Data  Sheet  and  distribution  list  only) 

Air  Force  Scientific  Adviser 
Director  Trials 

Aeronautical  and  Maritime  Research  Laboratory 

Director 

Chief  of  Airframes  and  Engines  Division 

Research  Leader  Structural  Mechanics 

Task  Manager  -  L.  Molent 

Author (s):  R.  Jones,  M.  Knopp,  J.  Price  (Monash) 

Others:  RLSI 

D.  Graham 
C.  Wang 
K.  Watters 
K.  Walker 
B.  Aktepe 
T.  Dickinson 

DSTO  Library 
Library  Fishermens  Bend 
Library  Maribymong 
Library  Salisbury  (2  copies) 

Australian  Archives 

Library,  MOD,  Pyrmont  (Doc  Data  sheet  only) 


L.F.R.  Rose 
A.  Mills 
W.  Hu 
W.  Guo 

CF IFOSTPTLO 


|  shared  copy 


Capability  Development  Division 

Director  General  Maritime  Development  (Doc  Data  Sheet  only 
Director  General  Land  Development  (Doc  Data  Sheet  only 
Director  General  C3I  Development  (Doc  Data  Sheet  only) 

Army 

ABCA  Office,  G-l-34,  Russell  Offices,  Canberra  (4  copies) 

SO  (Science),  DJFHQ(L),  MILPO  Enoggera,  Queensland  4051  (Doc  Data  Sheet 
only) 

NAPOC  QWG  Engineer  NBCD  c/-  DENGRS-A,  HQ  Engineer  Centre  Liverpool 
Military  Area,  NSW  2174  (Doc  Data  Sheet  only) 

Air  Force 

ASI1-DGTA  (2  copies) 

RAAF  TLO IFOSTP 

Intelligence  Program 

DGSTA  Defence  Intelligence  Organisation 

Corporate  Support  Program  (libraries) 

OIC  TRS,  Defence  Regional  Library,  Canberra 

Officer  in  Charge,  Document  Exchange  Centre  (DEC),  (Doc  Data  Sheet  and 
distribution  list  only) 

*US  Defence  Technical  Information  Center,  2  copies 
*UK  Defence  Research  Information  Centre,  2  copies 
*Canada  Defence  Scientific  Information  Service,  1  copy 
*NZ  Defence  Information  Centre,  1  copy 
National  Library  of  Australia,  1  copy 

UNIVERSITIES  AND  COLLEGES 

Australian  Defence  Force  Academy 
Library 

Head  of  Aerospace  and  Mechanical  Engineering 
Senior  Librarian,  Hargrave  Library, 

Monash  University  Librarian, 

Flinders  University 

OTHER  ORGANISATIONS 


NASA  (Canberra) 
AGPS 


OUTSIDE  AUSTRALIA 


ABSTRACTING  AND  INFORMATION  ORGANISATIONS 

INSPEC:  Acquisitions  Section  Institution  of  Electrical  Engineers 
Library,  Chemical  Abstracts  Reference  Service 
Engineering  Societies  Library,  US 

Materials  Information,  Cambridge  Scientific  Abstracts,  US 
Documents  Librarian,  The  Center  for  Research  Libraries,  US 

INFORMATION  EXCHANGE  AGREEMENT  PARTNERS 

Acquisitions  Unit,  Science  Reference  and  Information  Service,  UK 
Library  -  Exchange  Desk,  National  Institute  of  Standards  and  Technology,  US 
National  Aerospace  Laboratory,  Japan 
National  Aerospace  Laboratory,  Netherlands 

OTHER  ORGANISATIONS 

NAVAIRSYSCOM,  Structures  Division,  PATUXENT  RIVER,  MD  -  Atten.  Mr  Don 
Polakovics 

Boeing,  McDonnell  Aircraft  and  Missile  Systems  Company,  St.  Louis  -  Atten.  Eric 
Meyer 

National  Research  Council,  Institute  for  Aerospace  Research, 

-  Library 

-  Atten.  D.  Simpson 

Canadian  Forces,  DTA  3-8 
Bombardier  Services  -  Atten.  Y.  Richards 


SPARES  (5  copies) 

Total  number  of  copies: 


74 


Page  classification:  UNCLASSIFIED 


DEFENCE  SCIENCE  AND  TECHNOLOGY  ORGANISATION 
DOCUMENT  CONTROL  DATA 


2.  TITLE 


Stress  and  Strain  Estimation  at  Notches  in  Aircraft  Structures 


4.  AUTHOR(S) 

R.  Jones,  M.  Knopp,  J.  Price  and  L.  Molent 


1.  PRIVACY  MARKING/CAVEAT  (OF 
DOCUMENT) 


3.  SECURITY  CLASSIFICATION  (FOR  UNCLASSIFIED  REPORTS 
THAT  ARE  LIMITED  RELEASE  USE  (L)  NEXT  TO  DOCUMENT 
CLASSIFICATION) 

Document  (U) 

Title  (U) 

Abstract  (U) 


5.  CORPORATE  AUTHOR 

Aeronautical  and  Maritime  Research  Laboratory 
PO  Box  4331 

Melbourne  Vic  3001  Australia 


6a.  DSTO  NUMBER 

DSTO-GD-0196 

COE-SM-98-Ol 


6b.  AR  NUMBER 
AR-01 0-652 


6c.  TYPE  OF  REPORT 
General  Document 


8.  FILE  NUMBER 

9.  TASK  NUMBER 

10.  TASK  SPONSOR 

11.  NO.  OF  PAGES 

Ml/9/418 

AIR  98/092 

DGTA 

30 

7.  DOCUMENT  DATE 
October  1998 


12.  NO.  OF 

REFERENCES 

33 


13.  DOWNGRADING/DELIMITING  INSTRUCTIONS 


15.  SECONDARY  RELEASE  STATEMENT  OF  THIS  DOCUMENT 


14.  RELEASE  AUTHORITY 


Chief,  Airframes  and  Engines  Division 


Approved  for  public  release 


OVERSEAS  ENQUIRIES  OUTSIDE  STATED  LIMITATIONS  SHOULD  BE  REFERRED  THROUGH  DOCUMENT  EXCHANGE  CENTRE,  DIS  NETWORK  OFFICE, 
DEPT  OF  DEFENCE,  CAMPBELL  PARK  OFFICES,  CANBERRA  ACT  2600 


16.  DELIBERATE  ANNOUNCEMENT 
No  Limitations 


17.  CASUAL  ANNOUNCEMENT  Yes 


18.  DEFTEST  DESCRIPTORS 

Fatigue,  Stress  Analysis,  Fracture,  Damage  Tolerance,  Inelastic  Stress,  Plasticity  Theory 


19.  ABSTRACT 

To  maintain  the  continued  airworthiness  of  military  aircraft  it  is  essential  that  the  fatigue  behaviour  of 
components  subjected  to  complex  multi-axial  stress  conditions  be  both  understood  and  predicted.  This 
topic  is  extremely  complex.  Numerous  fatigue  failure  criteria  ranging  from  the  purely  empirical  to  the 
theoretical  have  been  proposed.  These  criteria  rely  on  the  estimation  of  stress  and  strain  at  fatigue  critical 
locations.  This  interim  report  focused  on  possible  approaches  which  may  be  applicable  to  both  low-  and 
high-cycle  fatigue  regimes.  It  discusses  the  relative  advantages  of  the  Neuber  and  the  Glinka  methods  for 
calculating  localised  notch  strains  as  compared  to  the  results  from  finite  element  analysis.  These  former 
techniques  are  at  the  core  of  several  sequence  accountable  crack  initiation  prediction  models,  some  of 
which  are  used  in  the  life  assessment  of  RAAF  aircraft.  Thus  the  accuracy  of  these  techniques  directly 
impact  the  estimated  fatigue  life  of  these  aircraft. 


Page  classification:  UNCLASSIFIED 


